Non-compactness and infinite number of conformal initial data sets in high dimensions

On any closed Riemannian manifold of dimension greater than $7$, we construct examples of background physical coefficients for which the Einstein-Lichnerowicz equation possesses a non-compact set of positive solutions. This yields in particular the existence of an infinite number of positive solutions in such cases